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Abstract 

We quantize a multidimensional SDE (in the Stratonovich sense) 
by solving the related system of ODE^s in which the d-dimensional 
Brownian motion has been replaced by the components of functional 
stationary quantizers. We make a connection with rough path theory 
to show that the solutions of the quantized solutions of the ODE con- 
verge toward the solution of the SDE. On our way to this result we 
provide convergence rates of optimal quantizations toward the Brow- 
nian motion for ^-Holder distance, q > 2, in LP{F). 



Key words: Functional quantization, Stochastic Differential Equations, 
Stratonovich stochastic integral, stationary quantizers, rough path theory, 
Ito map. Holder semi-norm, p-variation. 



1 Introduction 

Quantization is a way to discretize the path space of a random phenomenon: 
a random vector in finite dimension, a stochastic process in infinite dimen- 
sion. Optimal Vector Quantization theory (finite-dimensional) random vec- 
tors finds its origin in the early 1950's in order to discretize some emitted 
signal (see [IQ]). It was further developed by specialists in Signal Processing 
and later in Information Theory. The infinite dimensional case started to 
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be extensively investigated in the early 2000's by several authors (see e.g. 

m, 0, m, ca, etc). 

In [20] , the functional quantization of a class of Brownian diffusions has 
been investigated from a constructive point of view. The main feature of 
this class of diffusions was that the diffusion coefficient was the inverse of the 
gradient of a diffeomorphism (both coefficients being smooth). This class 
contains most (non degenerate) scalar diffusions. Starting from a sequence of 
rate optimal quantizers, some sequences of quantizers of the Brownian diffu- 
sion are produced as solutions of (non coupled) ODE^s. This approach relied 
on the Lamperti transform and was closely related to the Doss-Sussman rep- 
resentation formula of the flow of a diffusion as a functional of the Brownian 
motion. In many situations these quantizers are rate optimal (or almost rate 
optimal) i.e. that they quantize the diffusion at the same rate 0((log A^)^2) 
as the Brownian motion itself where N denotes the generic size of the quan- 
tizer. In a companion paper (see [27]), some cubature formulas based on 
some of these quantizers were implemented, namely those obtained from 
some optimal product quantizers based on the Karhunen-Loeve expansion 
of the Brownian motion, to price some Asian options in a Heston stochastic 
volatility model. Rather unexpectedly in view of the theoretical rate of con- 
vergence, the numerical experiments provided quite good numerical results 
for some "small" sizes of quantizers. Note however that these numerical 
implementations included some further speeding up procedures combining 
the stationarity of the quantizers and the Romberg extrapolation leading to 

3 

a 0((log A^)^2) rate. Although this result relies on some still pending con- 
jectures about the asymptotics of bilinear functional of the quantizers, it 
strongly pleads in favour of the construction of such stationary (rate optimal) 
quantizers, at least when one has in mind possible numerical applications. 

Recently a sharp quantization rate {i.e. including an explicit constant) 
has been established for a class of not too degenerate 1-dimensional Brow- 
nian diffusions. However the approach is not constructive (see ^). On 
the other hand, the standard rate 0((logA^)~2) has been extended in [22| 
to general d-dimensional Ito processes, so including d-dimensional Brown- 
ian diffusions regardless of their ellipticity properties. This latter approach, 
based an expansion in the Haar basis, is constructive, but the resulting 
quantizers are no longer stationary. 

Our aim in this paper is to extend the constructive natural approach ini- 
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tiated in |20) to general d-dimensional diffusions in order to produce some 
rate optimal stationary quantizers of these processes. To this end, we will 
call upon some seminal results from rough path theory, namely the continu- 
ity of the Ito map, to replace the "Doss-Sussman setting" . In fact we will 
show that if one replaces in an SDE (written in the Stratonovich sense) 
the Brownian motion by some elementary quantizers, the solutions of the 
resulting ODE^s make up some rough paths which converge (in p-variation 
and in the Holder metric) to the solution of the SDE. We use her the rough 
path theory as a tool and we do not aim at providing new insights on this 
theory. We can only mention that these rate optimal stationary quantizers 
can be seen as a new example of rough paths, somewhat less "stuck" to a 
true path of the underlying process. 

This work is devoted to Brownian diffusions which is naturally the promi- 
nent example in view of applications, but it seems clear that this could be 
extended to SDE driven e.g. by fractional Brownian motions (however our 
approach requires to have an explicit form for the Karhunen-Loeve basis as 
far as numerical implementation is concerned). 

Now let us be more precise. We consider a diffusion process 

dXt = b{t,Xt)dt + a{t,Xt)odWt, Xo = xGM'^, [0,r], 

in the Stratonovich sense where h : [0, T] x M'^ — )- M"^ and a : [0, T] x 
^ A4{dxd) are continuously differentiable with linear growth (uniformly 
with respect to t) and W = (VFj)tg[o,T] is a d-dimensional Brownian motion 
defined on a filtered probability space {Q,A,F). (The fact that the state 
space and W have the same dimension is in no case a restriction since our 
result has nothing to do with ellipticity) . 

Such an SDE admits a unique strong solution denoted = (Xf )^g[o,T] 
(the dependency in ^ will be dropped from now to alleviate notations). The 
M"^- valued process X is pathwise continuous and suptg[o^T] l^tl ^ U'{¥), r > 
(where | . | denotes the canonical Euclidean norm on M'^). In particular X 
is bi-measurable and can be seen as an L''(P)-Radon random variable taking 
values in the Banach spaces (L^ jgd) I • II^ ) where = L^a{[0, T],dt) and 

Iffli^ — ( So \9i't)\^dt^ ^ denotes the usual L^-norm when p£ [1, oo). 

For every integer > 1 , we can investigate for X the level N [U (P) , )- 
quantization problem for this process X, namely solving the minimization 
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of the L''(P)-mean L^j^^-quantization error 

eN,riX, L''') := mm^e]\f,r{ct, X, L''^), a C L^j^d, carda < (1-1) 

where eN,rioi, X, L^) denotes the L'^-mean quantization error induced by a, 
namely 

e^rid, X, V) := ( IE min \ X — a 

The use of "min" in (jl.ip is justified by the existence of an optimal 
quantizer solution to that problem as shown in [3l [13] in this infinite dimen- 
sional setting. The Voronoi diagram associated to a quantizer a is a Borel 
partition {Ca{a))a£a such that 

Caia) C \x£ LZ^^Wx - al^p <min|3;-6|^P [ 

and a functional quantization of X by a is defined by the nearest neighbour 
projection of X onto a related to the Voronoi diagram 

— I^al{XeC4«)}- 

In finite dimension (when considering M'^-valued random vectors instead 
of L^jg^-valued processes) the answer is provided by the so-called Zador 
Theorem which says (see [I0|) that if E\X\''+^ < +00 for some 5 > and if 
g denotes the absolutely continuous part of its distribution then 

N^eN,r{X,R'^) ^ Jr,d\\g\\K_ as N ^ 00 (1.2) 

d+r 

where J^^d is finite positive real constant obtained as the limit of the nor- 
malized quantization error when X = U{[0, 1]). This constant is unknown 
except when d = 1 or d = 2. 

A non- asymptotic version of Zador 's Theorem can be found e.g. in |22j : 
for every r, 6 > there exists a universal constant Cr^s > and an integer 
Nr^S ^ such that, for every random vector ^,A,¥) — )• M"', 

ViV > Nr,S, eNAX,^'^) < Cr,s\\X\\r+sN--^. 



min \X — a\ip 
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The asymptotic behaviour of the L'*(P)-quantization error of sequences of 
L'^-optimal quantizers of a random vector X when s > r has been extensively 
investigated in [13] and wih be one crucial tool to establish our mains results. 

In infinite dimension, the case of Gaussian processes was the first to 
have been extensively investigated, first in the purely quadratic case (r = 
p = 2): sharp rates have been established for a wide family of Gaussian 
processes including the Brownian motion, the fractional Brownian motions 
(see |18[ I19j). For these two processes sharp rates are also known for p G 
[1,00] and rG (0, 00) (see [1]). More recently, a connection between mean 
regularity of t Xt (from [0,T] into L^(P)) and the quantization rate has 
been established (see |22j): if the above mapping is /i-H61der for an index 
IJ,€ (0, 1], then 



Based on this result, some universal quantization rates have been obtained 
for general Levy processes with or without Brownian component some of 
them turning out to be optimal, once compared with the lower bound esti- 
mates derived from small deviation theory (see e.g. or [5]). One impor- 
tant feature of interest of the purely quadratic case is that it is possible to 
construct from the Karhunen-Loeve expansion of the process two families of 
rate optimal (stationary) quantizers, relying on 

- sequences {a^^'^^°'^^)N>i of optimal product quantizers which are rate 
optimaH.e. such that eAr^r'(«''^'*5 -'^i L'^) = 0{eN,2iX, L"^)) (although not with 
a sharp optimal rate). 

- sequences of true optimal quantizers (or at least some good nu- 
merical approximations) (a*-^'*^)7v>i i.e. such that eN,r{oi^^'*\X,Lr') = 



We refer to Section [2.11 below for further insight on these objects (both 
being available on the website www. quantize .math-fi . com). 

The main objective of this paper is the following: let (a^)Ar>i denote a 
sequence of rate optimal stationary (see (12.81) further on) quadratic quan- 
tizers of a (i'-dimensional standard Brownian motion W = iW^ , . . . ^W^). 
Define the sequence = {x^)n=i,...,N, N > 1, oi solutions of the ODE's 



e;v,r(X,iP) = 0((logiV)-^) 



pG (0,r). 



eiv,2(^,i'). 



x^{t) = x+ b{x^{s))ds+ a{xi:{s))da^{s), n=l,... 





N. 
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Then, the finitely valued-process defined by 



N 
n=l 

converges toward the diffusion X on [0, T] (at least in probability) as N ^ 
oo. This convergence will hold with respect to distance introduced in the 
rough path theory (see [251 El El El |26] ) which always implies convergence 
with respect to the sup norm. The reason is that our result will appear as 
an application of (variants of the) the celebrated Universal Limit Theorem 
originally established by T. Lyons in |25) . The distances of interest in rough 
path theory are related to the |-Holder semi-norm or the g-variation semi- 
norm both when q> 2 defined for every xG C([0, TjjR'^) by 

II II rr.- l^i't) ~ 

\\x\\q,Hol = Tp sup J < +00, 

0<s<t<T |t-s|9 

and 

^^^q,[0,T]{x) ■■= SUp|(^ ^ \xiti+i)-x{ti)\'^y 

o<e<k~i 

,0 <to <ti < <tk <T,k>l^ < +00 

respectively. Note that 

- a:;(0)||sup < Varpjo.T] (a;) < ||x||p,_f/oZ- 

From a technical viewpoint we aim at applying some continuity results 
established on the Ito map by several authors (see e.g. [HI [251 El [E]) that 
is the continuity of a solution x of the ODE (in a rough path sense) 

dxt = f{xt)dyt, xo = x{0), 

as a functional of y. However, the above (semi-)norms associated to a func- 
tion X are not sufficient and the natural space to define such rough ODE 
is not the "naive" space of paths but a space of enhanced paths, which in- 
volves in the case of a multi-dimensional Brownian motion the mutual Levy 
areas of its components. Convergence in this space is defined by considering 
appropriate |-Holder and p-variation semi-norms to both the d-dimensional 
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path and the related (pseudo-)Levy areas (wit different values of q and p, 
see Section [3]). Our application to quantized SDE^s will make extensively 
use the fact that our functional quantizations of the Brownian motion W 
will all satisfy a stationary assumption i.e. 

W = E{W\a{W)) 

so that we will extend the Kolmogorov criterion satisfied by W to its func- 
tional quantizers W for free. This approach is rather straightforward and 
its field of application seems more general than our functional quantization 
purpose: thus the piecewise affine interpolations of the Brownian motion 
obviously satisfy such a property (see Appendix). 

The paper is organized as follows. In Section 2 we provide some short 
background on functional quantization as well as preliminary elementary re- 
sults on stochastic integration with respect to a stationary functional quan- 
tizer of a d-dimensional standard Brownian motion. In Section 3, we define 
a quantized approximation scheme of an SDE (in the Stratonovich sense) 
driven by a standard Brownian motion by its functionally quantized coun- 
terpart which turns out to be a system of (non-coupled) ODE^s. To this 
end we recall some basic facts on rough path theory, in particular the notion 
of convergence we need to define on the so-called multiplicative functionals 
involved in the continuity of the Ito map which, when dealing with Brown- 
ian motion amounts, to some convergence in Holder semi- norm of the naive 
path as well as, roughly speaking, the running (pseudo-)Levy areas of its 
components. In Section 4 and 5, we establish successively the convergence 
in the Holder distance of sequences of optimal stationary quantizations W 
of the Brownian motion toward W: Section 4 is devoted to the conver- 
gence of the "regular" paths whereas Section 5 deals with the convergence 
of the running (pseudo-)Levy areas (and to the global convergence of the 
couple). In both cases we provide some convergence rate in the (logA^)"", 
aG (0, ^) scale which is the natural scale for such convergences since opti- 
mal functional quantizations of the Brownian motion are known to converge 
at a (log A^)~ 2-rate for most usual norms (like quadratic pathwise norm on 
LH[0,T],dt)). 

Notations: • For every d> I, one denotes ^ = (^^, . . . ,^'^) a row vector of 
M'^. A4{d X d) will denote the set of square matrices with d lines. 
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• I . I denotes the canonical Euclidean norm on W^. 

• We denote {T^)t>o the augmented natural filtration of a process X = 

{Xt)t>o (so that it satisfies the usual conditions). 

• For a bounded function / : [0,T] R'^, \\ f ||sup := sup^gp -r] \ f{t)\. If / is 

a Borel function and pe [l,+oo), := ( fl' \f(t)\PdtY . 

• For an M*^- valued bi-measurable process X and pG [l,+oo), we denote 

• We denote = A; = 0, . . . , 2", the uniform mesh of the interval [0, T], 
r>Oand/^ = fe = 0,...,2--l. 

• [x\ denotes the lower integral part of x € M. 

• Let (a„)„>o and (6„)„>o be two sequences of real numbers: a„ ~ 6„ if 
an = bn + o{bn) and a„ x 6„ if a„ = 0(6„) and 5„ = 0(a„). 

2 Background and preliminary results on func- 
tional quantization 

2.1 Some background on functional quantization 

Functional quantization of stochastic processes can be seen as a discretiza- 
tion of the path-space of a process and the approximation (or coding) of a 
process by finitely many deterministic functions from its path-space. In a 
Hilbert space setting this reads as follows. 

Let {H, (•,•)) be a separable Hilbert space with norm | • | and let X : 
{Q,A,f) H he a random vector taking its values in H with distribution 
Px- Assume the integrability condition 

E\X\'^<+oo. (2.3) 

For N>1, the L^-optimal iV-quantization problem for X consists in mini- 
mizing 

min \ X — a 

over all subsets a C H with card(a) < A^. Such a set a is called A^- 
codebook or A^-quantizer. The minimal quantization error of X at level A^ 
is then defined by 

{X, H) := inf | (E min \X - a\^f'^ : a ^ H, card(a) < n\ . (2.4) 



Emin X — a 
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For a given A^-quantizer a one defines an associated nearest neighbour pro- 
jection 

VTa := ^alc,(a) 



and the induced a- (Voronoi) quantization of X by setting 

X- :=7r„(X), 



(2.5) 



where {Ca(a) : a G a} is a Voronoi partition induced by a, that is a Borel 
partition of H satisfying 



Ca{oL) C {x G -ff : |x — a| = min |x — 6|} 



(2.6) 



for every a G a. Then one easily checks that, for any random vector X : 
a C H, 

E|X > E\X - = E min|X - 



SO that finally 

en{X,H) = inf<! |X-(7(X)| 
= inf 



q : H ^-^'^ card(g(F)) < N 



X - Y 



LH 



Y:{Vt,A) "4- H,caxd{Y{VL)) < n\.{2.7) 



A typical setting for functional quantization is H = := ^^([0, l],dt) 
(equipped with {f,g)2 ■= fg{t)dt and 1/1^2^ := y/{fJh). Thus any (bi- 
measurable, real-valued) process X = iXt)t<=[o,T] defined on a probability 
space {Q,,A,F) such that 



/ E{X'^)dt < +00 

^0 



is a random variable X : {Q,A,] 



L?. But this Hilbert setting is not 



the only possible one for functional quantization (see e.g. [21], [12], |5], etc) 
since natural Banach spaces like L^([0, T], (ii) or C([0,T],M) are natural 
path-spaces. 

In the purely Hilbert setting the existence of (at least) one optimal N- 
quantizer for every integer > 1 is established so that the infimum in (12. 4p 
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holds as a minimum. A typical feature of this quadratic Hilbert framework 
is the so-called stationarity (or self-consistency) property satisfied by such 
an optimal A^-quantizer a^^'*^: 

1"<"'*'=E(X|X-^"'*^ (2.8) 

This property, known as stationarity, will be used extensively throughout 
the paper. 

This existence property holds true in any reflexive Banach space and 
path spaces (see [12] for details). 



2.2 Constructive aspects of functional quantization of the 
Brownian motion 

2.2.1 Karhunen-Loeve basis (d = 1) 

First we consider a scalar Brownian motion (Wf)tg[o,T] on a probability space 
(ri,^,P). The two main classes of rate optimal quantizers of the Brownian 
motion are the product optimal quantizers and the true optimal quantizers. 
Both are based on the Karhunen-Loeve expansion of the Brownian motion 
given by 

Wt = Y,Vh^ke'^{t) (2.9) 



fc>i 



where, for every k > 1, 



and 



[ Wtsm{t/y% 
Jo 



The sequence (e^)fc>i is an orthonormal basis of L^. The system (A^, e^)fc>i 
can be characterized as the eigensystem of the symmetric positive trace class 
covariance operator of / H> (t H> J^{sAt) f{s)ds) = (t i-)- E(</ | W>2 Wt). 
In particular this implies that the Gaussian sequence (^a;)a:>i is pairwise un- 
correlated hence i.i.d., M{0; l)-distributed. The Karhunen-Loeve expansion 
of W plays the role of PC A of the process: it is the fastest way to exhaust 
the variance of W among all expansions on an orthonormal basis. 
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The convergence of the series in the right hand side of (j2.9p holds in 
for every cj G and ¥(dLjj)-a.s. for every tG [0,T]. In fact this convergence 
also holds in L^(P) and F{duj)-a.s. for the sup norm over [0,T]. The first 
convergence follows from Theorem I4.3f a) further on applied with X = W 
and Qn = cr(?i, • • • > ^jv) and the second one follows e.g. from [21j ¥{duj)-a.s.. 
In particular the convergence holds in Lp'[dF(^dt) or equivalently in L?2 (IP)- 

T 

Note that this basis has already been used in the framework of rough path 
theory for Gaussian processes, see e.g. [2l[7l[8]. 

2.2.2 Optimal product quantization {d > 1) 

> The one- dimensional case d = 1. The previous expansion of the Brownian 
motion suggests to define a product quantization of W at level N by 



{m,...,N^) 



where A^i , . . . , A'^^ are non zero integers satisfying A'^i • • • A'^^, < and ^{^^ , • • • , Cl^^ 
are optimal quadratic quantizations of .^i, . . . , . The resulting (squared) 
quadratic quantization error reads 

\\W_Wi^u-,N,)^^2 (2.12) 

k>l « 

An optimal product N -quantization w'^'^^"^ is obtained as a solution to 
the following integral bit allocation optimization problem for the sequence 
(A^fc)fc>i: 



mm 



{||Vr-I?(^i'-'^i)||2, A^i,...,A^^ > 1, A^r--A^L < A^, L>l} (2.13) 

(see [18] for further details and [27] for the numerical aspects). It is es- 
tablished in |18j (as a special case of a more general result on Gaussian 
processes) that 

L\\w - w^'P''°% X (log A^)-i (2.14) 

Furthermore, the critical dimension L = Lw[N) satisfies Lw[N) ~ log A''. 
Numerical experiments carried out in [27J show that 
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with ^ 0.5 (at least up to iV < 10000). 

It is possible to get a closed form for the underlying optimal product 
quantizers . First, note that the normal distribution on the real line 
being log-concave, there is exactly one stationary quadratic quantizer of full 
size M for every M > 1 (hence it is the optimal one). So, let > 1 and let 
(-^fc)fc>i denote its optimal integral bit allocation for the Brownian motion 
W. For every Nj, > 1, we denote by /S^^*) := {/3|^*\ 1 < ik < Nj,} the 
unique optimal quantizer of the normal distribution: thus a(0) = {0} by 
symmetry of the normal distribution. Then, the optimal quadratic product 
iV-quantizer a^-^™"' (of "true size" A'^i X • • • X -/V^^j^^jv) — ^) can be described 
using a multi-indexation as follows: 

<rl,...)(*) = E C''^ v^^rw, n,e{i,..., N,}, k > 1. 

k>l 

These sums are in fact all finite so that all the functions a^.'^^°'^ are C°° 
with finite variation on every interval of R+ . 

Explicit optimal integral hit allocations as well as optimal quadratic 
quantizations (quantizers and their weights) of the scalar normal distribution 
are available on the website [28]. Note for practical applications that this 
optimal product quantization is based on 1-dimensional quantizations of 
small size of the scalar normal distribution AA(0; 1). This kind of functional 
quantization has been applied in [27J to price Asian options in a Heston 
stochastic volatility model. 

> The d-dimensional case. Assume now W = {W^, . . . , W^) is a d-dimensional 
Brownian motion. Its optimal product quantization at level > 1 will be 
defined as the optimal product quantization at level [A^^J of each of its d 
components. 

[> Additional results on optimal vector quantization of the normal distri- 
bution on M.'^. We will extensively make use of the distortion mismatch 
result established in [13] that we recall here only in the d-dimensional Gaus- 
sian case. Let Z be an AA(0; 1^) random vector and let be an optimal 
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quadratic quantizer at level N of Z (hence of size N). Then 

{i) ypG{0,2 + d),yN>l,\\Z-Z'^^\\p<Cz,pN-2, (2.15) 

(ii) VpG [2 + (i,+oo), Vr?G (0,d + 2),ViV > 1, 

WZ-Z'^^Wp < Cz,p,^N-'-^ (2.16) 

where Cz,p and Cz,p,r] are two positive real constants. 
2.2.3 Optimal quantization (d = 1) 

It is established in |18] (Theorem 3.2) that the quadratic optimal quantiza- 
tion of the one-dimensional Brownian motion reads 




(2.17) 



where, for every integer d > 1, (d = Proj^^(Ty) ~ AA(0; Diag(Ai, . . . , A^)) 
with Ed := M-span {sin (./"v/Ai) , • • • , sin (./V^) } and is an optimal 
quadratic quantization of (d at level (or of size) N. 

If one considers an optimal quadratic A^-quantizer = {/3^, n = 
1, . . . , iV} C M'^M' W of the distribution M (O; Diag(Ai, . . . , Xdw{N))) (« Pri- 
ori not unique) 

dw{N) 

«^ni)= 5] (/3f ))'^V^er(t), n = l,...,N. 

k=l 

Once again this defines a C°° function with finite variation on every interval 
of M+. 

A sharp rate has been obtained in [19] for the resulting optimal quanti- 
zation error 

\\W -W^'°P%^Tc°P\logN)-^2 as ^ oo (2.18) 
where c°p^ = ^ 0.4502. 

W TT 

The true value of the critical dimension dw{]^) is unknown. A conjecture 
supported by numerical evidences is that dw{]^) ~ log A^. Recently a first 
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step to this conjecture has been estabhshed in [23] by showing that 

N log(iV) - 2 

Large scale computations of optimal quadratic quantizers of the Brown- 
ian motion have been carried out (up to = 10 000 and d = 10). They are 
available on the website [28] . 

In the d-dimensional setting, several definitions of an optimal quantiza- 
tion of the Brownian motion W = {W^ , . . . jW^) can be given. For our 
purpose, it is convenient to adopt the following one: 



l<j<d 

Its property of interest is that this definition preserves the componentwise 
independence as well as a stationarity property (see below) since 

1 1 

e(w' I W^'°P') =e(w'\ Wi^^^^'°^'^ = wi^^''^'"^'^ i = l,...,d. 

2.2.4 Wiener like integral with respect to a stationary functional 
quantization (d = 1) 

Both types of quantizations defined above share an important property of 
quantizers: stationarity. 

Definition 2.1. Let a G L"^, a ^ $, be a quantizer. The quantizer a 
is stationary for the (one- dimensional) Brownian motion W if there is a 
Voronoi quantization W := W"' induced by a such that 

W = E{W\a{W)) a.s. (2.19) 

where E( . |^) denotes the functional conditional expectation given the a -field 
Q on L?2 (IP) (see Appendix) and aiW) is the a-field spanned by W. 

T 

Note that if a is stationary for one Brownian motion, so it is for any 
Brownian motion since this stationarity property only depends on the Wiener 
distribution. 

In the case of product quantization W^'^^°^, this follows from the sta- 
tionarity property of the optimal quadratic quantization of the marginals 
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(see [18] or [27]). In the case of optimal quadratic quantization W^'"^^ this 
fohows from the optimahty of the quantization of Cdw{N) itself. 

We will now define a kind of Wiener integral with respect to such a 
stationary quantization 1^ of a one- dimensional W. So we assume that 
d = 1 until the end of this Section. 

First, we must have in mind that if W is an (J^j)-Brownian motion where 
the filtration {J^t)t>o satisfies the usual conditions, one can define the Wiener 
stochastic integral (on [0,r]) of any process ip £ L'^{[0,T] x i},B{[0,T]) (g) 
J-o,dt dP) with respect to W. The non-trivial case is when J-^ ^ Tt, 
typically when Tt = V ^ t G [0,^] where B and W are independent. 
One can see it as a special case of Ito stochastic integral or as an extended 
Wiener integral: if (<^(t, j)£Qx[o,t] denotes an elementary process of the 
form 

n 

•^{t.iS) := ^¥'A:(w)lsj^<t<Sft+i, = So < Si < • • • < Sn-1 < S„ = T 
fc=l 

where the random variables are J-"o-measurable (hence independent of 
W). Set 

n 
k=l 

Then, is an isometry from L?2 (P) into L^(J-"j,,P). Furthermore, one 

T 

easily checks that 

E <^(s, .)dWs I J-f ) = j\ {ip{s, .) I J-^) dWs 

where J-^ denotes the augmented filtration of W at time T. We follow 
the same lines to define the stochastic integral with respect to a stationary 
quantizer. Set for the same elementary process if 

n 

k=l 
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so that 



k=l 
n 



k=l 



E 



T 



ip{t,.)dWt\Toya{W) 



where we used that the (T-fields a{W) and J-q are independent since is a 
Borel function of W. As a consequence, 

WiMWl < WiMWl = \\ M^^e 

Hence, the hnear transformation extends into a hncar continuous map- 
ping on the whole set L^a (J'o,^)- Furthermore, one checks, first on el- 
ementary processes, then on {J^o,^) by continuity of the (functional) 

T 

conditional expectation, that 

We will denote from now on I^{ip){Lo) as an integral, namely 

Jo 

Now set as usual, for every te [0,r], 

/ (p{s,uj)dWs{Lo) := / lio^t]{s)(pis,uj)dWs{io). 
Jo Jo 

One checks using Jensen and Doob Inequality that. 



E sup 
te[o,T] 



< E sup 

te[o,T] 

< 4E 

'0 



ip{s, .)dWs 

/ ip'^{s,.)ds. 
Jo 



(2.20) 
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Furthermore, as soon as the underlying stationary quantizer a (such that 
W = T^") is made up with pathwise continuous elements, for every elemen- 
tary process its integral process 

t ri 

ip[s,.)dWs=Y,UWs,^,f^t-Ws,M) 
k=l 







pathwise continuous as well since W is a-valued. One classically derives, by 
combining this result with (j2.2Up and the everywhere density of elementary 
processes, that, for every if£ L^2 {J^o,^)i the process 

T 

ip{s, .)dWs 1 admits a continuous modification. 

J te[o,T] 

This is always this modification that will be considered from now on. As a 
matter of fact, if ifn denotes a sequence of elementary processes in (J^o,^) 

T 

converging to ip, i.e. satisfying 



E / {ip — ip-n)^ {s , .)ds — )• as n — )• oo. 
Jo 

It follows from (j2.20p that the convergence also holds in L r oo (^o i IP) • In 

T 

particular, there is a subsequence that converges T-a.s. for the || . ||sup which 

implies the existence of a continuous modification for / ip){s,ijj)dWs{u^). 

Jo 

Finally, using the characterization of functional conditional expectation 
(see Appendix), it follows that 

E (^^ vis, .)dWs , I Jo V a{W)^ = ^ vis, .)dWs. (2.21) 

Proposition 2.1. Let W be a (real-valued) Tt- standard Brownian motion. 
(a) For every ipe (J^o,^) 

T 

[ v{s,.)dWs = xl^y^ik I v{s,-)cos{s/^/X~k)ds (2.22) 
Jo V 1 Jo 



where '■= {W\eY)2/ ^'^^ independent, M{0;1)- distributed (see 
and 1^2. 10\) ) and independent ofip- 
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(b) Let W be a stationary quantization ofW. For every ip^ L?2 (-^o 

T 



{W\e 



k )2 



k>l V 



ip{s, .) cos{s/ \/ Xn)ds. (2.23) 



In particular if is a product quantization, then 



Proof, (a) Set for every ip€ {Fq,W), 



xl^yZCk^/h [ ^{s,.)dsm{s/y^k) (2.24) 
^ ^ k>i 

^ ^ k>l 



This defines clearly an isometry from L?2 {^q,^) into the Gaussian space 

T 

spanned by (Cn)n>i since 

10 



nJAv?) = |^E(ei)E (/5(s)-^cos(s/ 



/VAfc)(is ) = E / 



The last equality uses that the sequence ^y^cos(7r(/c — ^)t/T)j ^ is an 

orthonormal basis of L^. Finally, note that for every t£ [0,T], Jj,(l[o_t]) = 

•\/^ Sa:>i V^Cfc sin(t/-v/Afc) = l^t- This proves that Jj, = Ij, i.e. is but the 
(extended) Wiener integral with respect to W . 

(6) This follows by taking the (functional) conditional expectation of (j2.22p . 

2.2.5 Application to multi-dimensional Brownian motions {d > 2) 

Now we apply the above result to a componentwise (stationary) functional 
quantization of a multi-dimensional standard Brownian motion. 
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Proposition 2.2. Let W =: {W^, ■ ■ ■ , W^) denote a d- dimensional standard 
Brownian motion and let W := (W^ , . . . jW^) be a pathwise continuous 
stationary quantization of W (no optimality is requested here). Then, P- 
a.s., for every i ^ j, i, j G {1, . . . , d}, for every s,t£ [0, T], < s < t, 

E V; - W:)dWi I a(W)^ = j\wi - WDdWi 

Proof. All the components of W being independent, it is clear one can le- 
place (t(VF) by (t(VF% VK-'). Then, the stochastic integral j^WldWi coincides 
with the (extended) Wiener integral defined with respect to the filtration 
Ql^ := cj{J^^^ (it is clear that is a ^^^-standard Brownian motion 

still by independence). The result is then a straightforward consequence 

Remark. The above result still holds if one considers an additional "0*'"' com- 
ponent =t to the Brownian motion and to its functional quantization 
by setting W° = t as well. 

3 Convergence of quantized SDE''s: a rough path 
approach 

3.1 Prom Ito to Stratonovich 

An SDE 

dXt = b{t,Xt)dt + a{t,Xt)dWt, XoS L^,(P) 

where b : [0,T] x R'^ ^ R'^ and a : [0,T] x R'^ ^ M{d x q) are smooth 
enough functions (e.g. continuously differentiable with bounded differen- 
tials) and W = {Wt)te[o^T] is a (/-dimensional Brownian motion. First note 
that without loss of generality one may assume that q = d hy increasing 
the dimension of W or adding some identically zero components to X (no 
ellipticity like assumption is needed here). This SDE can be written in the 
Stratonovich sense as follows 

dXt = f{Xt)odWt, XoGL^,(P), (3.25) 

where, for notational convenience W = (W^, VF^, . . . , W^) stands for {t, Wt), 
Xt = {X^,X^, Xf) stands for {t, Xt) and / : [0, T]xR'^ ^ M{{d + 1) x 
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{d + 1)) (with fit, x) = {1,0,..., 0) as O*'^ line) is a differentiable function 
with bounded differentials. 

Following rough paths theory initiated by T. Lyons ([25]) and developed 
with many co-authors (see e.g. \26 \ \14: \ \T6 [ \26 \ \9\ for an introduction), one can 
also solve this equation in the sense of rough paths with finite p-variation, 
p > 2, since we know {e.g. from the former Kolmogorov criterion) that W 
a.s. does have finite ^-Holder norm, for any q > 2. Namely this means 
solving an equation formally reading 

dxt = f{xt)dyt, xoGM'^. (3.26) 

In this equation y does not represent the path (null at 0) yt = Wt{uj), 
t £ [0, T] itself but an enhanced path embedded in a larger space, also 
called geometric multiplicative functional lying on y with controlled |-Holder 
semi-norm, namely a couple y = {{yl^t)o<s<t<T, {ys,t)o<s<t<T) where ^ = 

can be identified with the path {yt) and 
{ylt)o<s<t<T satisfies, y^ ^ G M^'^+i)' for every < s < n < t < T and the 
following tensor multiplicative property 

ylt = ylu + yl,t + yl,u ® yi,f 

Different choices for this functional are possible, leading to different so- 
lutions to the above Equation (|3.26p . The choice that makes coincide a.s. 
the solution of (j3.25p and the pathwise solutions of (j3.26p is given by 

ylt = Wt{u) - Ws{u), ylt := { f {W^ - Wl) o dWi] {uj) (3.27) 

\J s / i,j=0,...,d 

SO that 

•J s,u ^ J u,t s,uJ' s,uj i j—Q^ ii 

The term y^^ is but the "running" Levy areas related to the components 
of the Brownian mtion W . The enhanced path of W will be denoted W 
(although we will keep the notation y in some proofs for notational conve- 
nience). One defines, for every (7 > 1, the |-Holder distance by setting 

Pqiy - x) = ||y^ - y^^\\q,Hol + IIy^ - X^||g/2,HoZ 

where ^ 

\\^^\\q/2,Hol ■=T^ sup 



0<s<t<T \t- sU 
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Remark. Likewise, when p G [2,3), one defines the p-variation distance 
between two such multiphcative functionals y, z is defined by 



^p{y,^) = Varpjo,T](y^ - z^) + Varp/2,[o,T](y^ - z^) 

where 



Vargjo,T](y^):=sup< 



^\yl,tej'] ,0<to<ti<---<tk <T,k>l\. 



The distance pq has been introduced in |24j although rough path theory 
was originally developed for the distance 5p in p-variation. Recently several 
authors came back to Holder distances pq (see e.g. [T6l l6t |9] ) . 

The following so-called universal limit theorem theorem (including vari- 
ants) describes the continuity of the so-called Ito map y i— t- x with respect 
to both 6p and pp-distances and will be the key for our main result. It was 
the starting point of rough path theory initiated by T. Lyons. Several state- 
ments (or improvements) can be found e.g. in [25 \ \14 \ \TE [ 1261 [9]. We state 
here some versions coming from [14] and |16j . 



Theorem 3.1. Let a£ (0, 1]. 

(a) (See JIB^) Let f : [0,T] xR'' ^ 7W((d + 1) x (d + 1), twice differentiahle 
with a hounded first differential and an a-Holder second differential. Suppose 
the multiplicative functional y satisfies Hy"*^ — x^Hg^/^o; + ||y^ — ^'^\\q/2,Hol < 
+00 for q£ (2, 2 + a). Then Equation \3.26\) has a unique solution starting 
at xq. 

When y = W(a;) (i.e. given by \3.21 )), the first component x.^ = x of 



the solution solution x = (x-'^,x^) a.s. coincides with (-^t('^))te[o,r]) solution 
to the SDE in the Stratonovich sense. 

Furthermore, the Ltd map y i— )• x zs continuous for the Holder pq distance 
(and locally Lipschitz in sense described in tiSf ). 

(6) (See 13 If f e C2([0,T] x R'^,M{{d + 1) x {d + 1)) is such that 
f.Vf is hounded with an a-Holder differential, then the conclusions of claim 
(a) still hold. 
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3.2 Quantization of the SDE and main result 

Let (a^)Ar>i denote a sequence of quantizers of the Brownian motion. Each 
is made up of N functions (or elementary quantizer) : [0, T] — )• M*^, 
n = 1,. . . ,N. For convenience a component "0" will be added accordingly 
to each elementary quantizer by setting an'^{t) = t (which exactly 
quantizes the function = t). We assume that every elementary quantizer 
is a continuous function with finite variation over [0,T]. The resulting 
Voronoi quantizer W = W^^ of W reads 

N 

Wt = ^a^(i)l{H/eC4a^)}' [0,r]. 

n=l 

Our aim is to approximate the diffusion process {xt)te[o,T] solution to 
the SDE (j3.25p by the solution of the equation 

dX^ = f{xf')dWu X^ = xo. 

as — )• oo. In fact, a less formal expression is available for the process X^ , 
namely 

N 
n=l 

where each is solution to the ODE 

dx^{t) = f{x^{t))da^{t), x^{0)=xo, n = l,...,N. (3.28) 

Note that X^ is a non- Voronoi quantization of (xt) (at level A^). The start- 
ing natural idea was to hope that X'^ converges to (xt) owing to the con- 
vergence of toward W. . . in an appropriate sense. Since we will use the 
above Theorem 13.11 we need to prove the convergence of the geometric func- 
tional related to W toward that of W. The quantity is formally 
defined by mimicking the definition of W, namely, for every {s,t) G [0,r], 
0<s<t<T, 

Wi'^(a;) :=Wt{u:) - Ws{u:), W^f (u;) := ( [\w: - WDdwA (u) 

\Js / i,j=0,...,d 

still with the convention wj^''^ = t. The integral must be understood in the 
usual Stieltjes sense. 
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Theorem 3.2. Let {W )n>i be a sequence of stationary quadratic func- 
tional quantizers of the Brownian motion converging to W in (P). 

T 

Let f be like in claims (a) or (b) in Theorem \3.1[ Consider for every 
N > 1, the solutions of the quantized ODE 

dX^ = f{X?)dWi', N>1. 

as defined by hS. 28\) . LefK. and'K^ denote the enhanced paths of X, solution 
to 113. 25\) . and X^ respectively. Then, for every (2, 2 + a), 

Furthermore if r > ^ then 

p,(xLe-^J,X)^0. 

In view of what precedes this result is, as announced, a straightforward 
corohary of the continuity of the Ito map estabhshed Theorem 13. H once the 
convergence /9g(W^,W) in probabihty is estabhshed for any q€z (2,3). A 
shghtly more derailed proof is proposed at the end of Section 5. 

In fact we will prove a much precise statement concerning the Brownian 
motion since we will establish for every q > 2 the convergence in every LJ'{¥), 
< p < 00, of yOg(W^, W) with an explicit L^(P)-rate of convergence in the 
scale (logA^)-^ Be (0,1). 

These rates can be transferred to the convergence of the quantized SDE, 
conditionally to some events on which the Ito map is itself Lipschitz con- 
tinuous for the distances pq. Several results of local Lipschitz continuity 
have been established recently, especially in j6], [9], [16], [T7], although not 
completely satisfactory from a practical point of view. So we decided not to 
reproduce (and take advantage of) them here. 

The proof is divided into two steps: the convergence for the Holder semi- 
norm) of the regular path component is established in Section [J] (in which 
more general processes are considered) and the convergence of approximate 
Levy areas in Section [5] (entirely devoted to the Brownian case for the sake 
of simplicity). 

Remarks. • There is a small abuse of notation in the above Theorem since 
X^ is not a Voronoi quantizer of X: this quantization of X is defined on 
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the Voronoi partition (for the L^j^^-norm) induced by the quantization of 
the Brownian motion W. 

• The same results holds for the Brownian bridge, the Ornstein-Uhlenbeck 
process and more generally for continuous Gaussian semi-martingales that 
satisfy the Kolmogorov criterion. 

4 Convergence of the paths of processes in Holder 
semi-norm 

4.1 A general setting including stationary functional quan- 
tization 

In this section we investigate the connections between the celebrated Kol- 
mogorov criterion and the tightness of some classes of sequences of processes 
for the topology of ^-Holder convergence. In fact this connection is somehow 
the first step of the rough path theory, but we will look at it in a slightly 
different way. Whatsoever this naive pathwise convergence is not sufficient 
to get the continuity of the Ito map in a Brownian framework and we will 
also have to deal for our purpose with the multiplicative functional (see 
Section [5]) . 

But at this stage we aim at showing that when a sequence {Y^)n>i 
satisfies some "stationarity property" with respect to a process Y, several 
properties of Y can be transferred to the Y^ . Indeed, the same phenomenon 
will occur for the multiplicative function (see the next section). 

If Y satisfies the Kolmogorov criterion and (^Ar)Ar>i denotes a sequence 
of sub-cr-fields of A, then a sequence of processes defined by 

Y^ :=E{Y\g^), N>1, 

where the conditional expectation is considered in the functional sense (see Ap- 
pendix) is (C-tight and) tight for a whole family of topologies induced by 
convergence in ^-Holder sense. 

Definition 4.2. Let p > 1, 9 > 0. A process Y = (^)fg[o,T] satisfies the 
Kolmogorov criterion (Kp^g) if there is a real constant C^°'>0 such that 

Vs, te[0,r], E\Yt-Ys\P <C^°^\t- s\^+^ and YqG LP{¥). 
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Theorem 4.3. Let Y := (lt)te[o,T] a pathwise continuous process defined 
on {Q,A,F) satisfying the Kolmogorov criterion {Kp^). Let {Gn)n>i be a 
sequence of sub-a-fields of A. For every N > 1 set 

Y^ ■.= E{Y\gN). 

For every N > 1, Y^ has a pathwise continuous version satisfying 

yte[0,T], Yf = E{Yt\gj a.s. 

Furthermore, if one of the following conditions is satisfied 

(a) Qn C Qn+i, 

(b) There exists an everywhere dense subset D C [0, T] such that 

ytG[0,T], Yf^Yt. 

(c) \Y^ — Y\lt^ for some r > 1, 
then 

Vg>^, yp€[l,qe), \\Y -Y^\Up+\\Y -Y^Wg^Hol-^O. 

The proof of the theorem is a variant of the proof of the Kolmogorov 
criterion for functional tightness of processes. It consists in a string of several 
lemmas. For the following classical lemma, we refer to |14| (where it is stated 
and proved for semi- norms in variation). 

Lemma 4.1. Let x, y e C{[0,T],W^) and let q>l. Then 

(a) ||x - 2;(0)||sup < \\x\\gMol- 

(b) \\x + y\\q^Hol < \\x\\p,Hol + \\y\\q,Hol ifq>'^, 

(c) For every q> q' >1, < (2||x||sup)''"^' Ikll^,^^^,. 

(d) Claims (a)-(6)-(c) remain true with the p-variation semi-norm Var^ [o,T] 
instead of the --Holder semi-norm. 

J q 

Lemma 4.2. Letp€z [l,oo). LfY satisfies the Kolmogorov criterion [Kp^g) 
then, for every N > \, the process Y^ defined by Y^ = E(lt | Q^) has a 
continuous modification which is —-Holder continuous for every 9' £ (0,0) 
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(i.e. \\y \\£j.^Hoi < +00 a.s.). Furthermore, the sequence {Y is C- 

tight and for every 9' G (0,^), there exists a random variable Zqi G L^i^) 
such that 

nduj)-a.s. \\Y{oj)\\^^Hoi<Ze' (4.29) 

and 

ViV > 1, W{uj)\\2^^Hoi < I Gn){u;). (4.30) 

In particular, the sequence of Holder semi-norms {\\Y^ \\£j.^Hoi)n>i is 
-uniformly integrable. 

Remark. As a by-product of the proof we also get that 

where Cx^p^efi' is a finite real constant that only depends upon p, T, 6 and 
9' (and not on Y or the ir-fields Gpf). 

Proof. First it follows form the Kolmogorov criterion that for every > 1, 
Y^ admits a continuous modification which is ^-Holder for every 9' G (0, 9). 
Moreover the sequence {Y^)]\f>i is C-tight since every Y^ satisfies the same 
Kolmogorov criterion (Kp^g) and Yj/^ = E(yo|^iv) is tight on R (see [l], [29] 
p. 26). Now, let s,t£ [0, T], let m, n > 1 be two fixed integers. First note 
that 



sup lit— < max \Y,n + m — Y.n + m\ (4.31) 

t^lnT'] ^t^T 0<A;<2"+'"-l k+i 

s,t&[0,T\,t<s<t+i^ m>0 - - 

and 

max \Y, P < > IK 

n + m 1 ,n-\-m 

k=0 



For every 6' £ (0,6), set 



Ze' := I f E sup \Yt-Ys\\. (4.32) 

,n>0 s,t(£[0,T],t<s<t+^ 
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Taking the L^-norm in (j4.3ip yields 

^ / n>0 s,te[0,T],t<s<t+^ 
e' 

(2 \ ~ g' 
— I >^2"~ >^ II max \Y,n + m — Y,n + m 

T ^ ^ 0<fc<2"+™-l *k+i *fe 

^ n>0 £>0 - - 

On the other hand, owing to the the Kolmogorov criterion (Kp^g), 

' < > EIK 

0<A;<2" + '" — 1 k + 1 ''k ^ — ' ''k + 1 '■k 

- - k=0 



E 



^ 2n+mfjKol2-(,n+m){l+e)rp-{l+e) 
_ (jKolrp-{l+e)2-(n+m)e 



Hence 



\n>OTn>0 J 

where the finite real constant Cp Tfifi' only depends on T, 6 and 9' . On 
the other hand, for every 5 G [0,T], there exists a integer > 1 such that 
2-{i+n,) < ^ij. < 2-"^. Hence, 

r'^'sup \Yt - Ys\P < 2(1+"^)''' T-^' X sup \Yt - Y^l^ < Z^,. 

s,t£lO,T],t<s<t+5 s,t&[0,T],t<s<t+^ 

Consequently, for every s, tG [0, T], and every cjG ri, 
\Yt{uj)-Y,{io)\<Ze'{uj)\t-s\l^ 

i.e. 

\\Y{i^)\\^^Hoi<Ze'{^). 
Finally, it follows from Jensen's Inequality that for every s, tE Qn [0,T], 

¥{dw)-a.s. l^t^H - Y,^{uj)\ < E{Ze' \ aiv)H|t - sf. 
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In particular this means that, for every p> 1 and every 6' G {0,9), 

F{duj)-a.s. \\Y^{uj)\\£^^Hoi < HZe' I Qn){oo) < +oo 
and satisfies the LP-uniform integrabihty assumption. ^ 

Proof of Theorem[431 The sequence {Y'^)n>i being C-tight on {C{[0,T],R'^), \\ . ||sup), 

so is the case of the the sequence {Y^ , Y)n>i on (C([0, T],M?'^), \\ . ||sup) since 

the product topology coincides with the uniform topology. Let Q = 

limjv IP(yjv' y) denote a weak functional limiting value of (y^,y)jv>i- If 

H = {E^,E'^) denotes the canonical process on (C([0, T], M^'^), || . ||sup), it is 

clear that 0=2 = Py . 

> Convergence of the sup-norm. Assume that (c) holds: the functional 
y I— )• — 2/^(01^5, is continuous on (C([0, T], M^°'), || . ||sup)) consequently, 

- H^l^r = Q-a.s. i.e. Q = P(y,y) so that {Y^,Y) ^^'^""^ (Y,Y) as 

AT IP 

— 7- oo which simply means that ||y^ 

sup — ^ 0. On the other hand, it 
follows from Lemma 14.21 that, for every A > 1, 

WY"" - y||Lp < Cp,T inz^e' I Gn) + Zi:) a.s. 

(for a given fixed 9' G (0,^)) which implies that (Hl"^ — y||sup)Af>i is uni- 
formly integrable. Finally, 

E||y^ -yll^up ^ as A^oo. 

Assume that (6) holds: it follows that, for every ti, . . . , G D, one has 
{YJ^ , . . . , Y-i^) (Yj^ , • • • , Ytj, ), which in turn implies that the convergence 
{Y,^ , . . . ,Yt^ ,Yt„ . . . ,Yt,) A {Yt„...,Yt„Yt„...,Yt,). This means that 
Q and F(y,y) have the same finite dimensional marginals i.e. Q = P(yy). 
One concludes like in (c). 

If (a) holds, for every t£ [0,T], Y/^ — )• Yt F-a.s., so that (6) is satisfied. 

> Convergence of the Holder semi-norm. Let q > 1. As concerns the 
convergence of the |-Holder semi-norm, one proceeds as follows. Let q' G 

and set 9' := ^€ {0,9). It follows from Lemma[m6)-(c) that 

\\Y - Y'^Wg^Hol < 2^-^ ||y - y^lllup^ X {\\Y\\g,^Hol + \\Y''\\g',Hol)'^ . 
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Now let Z := Zq, be defined by gSS]). Then, 

\\y%',Hoi^\\y''\\q>,Hoi < z + (E(z|g^)). 

Hence, the sequence (||^||g',_H'o« + ||^^||g',_ffo«)Af>i> is tight since it is L^- 
bounded. On the other hand, ||y — I'^Hsup so that 
\\Y -Y^'Wq^Hoi^^^N 

Now let = I G (0,^). The same argument as above shows that ||y — 
y^\\q,Hoi < Z + E(Z|g^) where Z = is still given by g32|). As a 
consequence, {\\Y — \Tq hoi)^>^ uniformly integrable since, for every 

> 1, Jensen's Inequality implies 

\\Y -Y'^\\l^^,<2^-\Z^ + nZ'\G^)) 

AT LP 

which finally implies that — Y \\q^Hol — ^ 0. ^ 

4.2 Application to stationary quantizations of Brownian mo- 
tion: convergence and rates 

Theorem 4.4. (a) Let {W^)i\f>i be a sequence of stationary quadratic func- 
tional quantizers of a standard d- dimensional Brownian motion W defined 
by 12.11\) or ( f^. j7| j converging to W in a (purely) quadratic sense, namely 
II \W — W^\i2 II2 — )• as N ^ 00. Then, for every q > 2, 

VpG(0,oo), \\W -W^Wg^Hoi ^0 as N ^ 00. 

(b) Let q > 2. If, for every N > 1, W'^ is an optimal product quantization 
at level N. Then, for every p£ (0,c«), 

\W-W''\U,Hoi =o((logiV)-i'^^°(^(^-|)4)+"), V< 

p 

The proof of this Theorem is a consequence of the above Theorem 14. 3[ 
So we need to get accurate estimates for the increments of the processes 
W — W'^ . This is the aim of the following lemma. 

Lemma 4.3. Let p G [2, +00). Let , N > 1, denote a sequence of 
optimal product quadratic quantizers. For every p G (0, |) and every e G 
(0,3), for every s, t£ [0,T], s < t. 



'a > 0. 



{Wt - Ws) - {Wf - T^f ) < Cp,p,TAe\t - s\%log iV)-^5 



(4.33) 
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In particular, ifp£ (2,3), then 



<Cp,p,T,d|t-sr(logiVr^2 



(4.34) 



Proof. We may assume without loss of generality that we deal with a 
one-dimensional Brownian motion W , quantized at level N' = [A^dJ since 
everything is done component by component. Set for every A; > 1, '■= S,k~ 
where A'^i , . . . , Nf^ , . . . denotes the optimal bit allocation of an optimal 
product quadratic quantization at level A'^'. Keep in mind that for every 
k > L^{N'), Nk = 1 and that of course A^i • • • iVL^(Ar') < N' . The random 
vectors {Ck)k>i are independent and centered. 

It follows from the K-L expansion of W and its product quantization 
that 

(m - Ws) - - Wn = ^MeVit) - efis)) . 

k>l 

Then, it follows from the B.D.G. Inequality for discrete time martingales 
that 



{Wt-Ws)-{wf -lyf ) 



}N' 



< a 



Y^XkUe"^ {t) - eY {s)f 



k>l 



t - s[ 



, fc>l 



since, for every A; > 1, 



sm 



t - s 



cos 



t - s 



<^\t-s\"^Xl 



The random variables being an optimal quadratic quantization of 
the one-dimensional normal distribution for every k £ {1,. . ., L^{N')}, it 
follows from (12.16P that, there exists for every e £ (0,3), a constant Hp^^ 
such that 

vm>i, 11611, = II? -erii,<«p 



1 
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where ^™ denotes the (unique) optimal quadratic quantization at level m of 
a normally distributed scalar random variable ^. As a consequence, 



< Cp,T,e\t - s\ 



1 



,fc>i N, " 



> Temporarily assume that [2, 3). One may choose e so that 1 A = 1. 
Now, keeping in mind that L' := L^{N') ~ log A^' and < c/c~^ for a real 
constant c > 0, one gets 



fc ^ fc=l k>L' 



\ k=l k J 



Now, following e.g. [18], we know that the optimal bit allocation yields 

k=l k 



SO that, finally 



N' 



N' 



<Cp,p,T\t-s\P (log 



> Assume now that p£ [3, +oo) and e G (0,3). Set p = > 1 and q its 
conjugate exponent. Then, Holder Inequality implies 



.=1 Ni 



< 



P / L' 



_ . 1 pp 



\k=l 



We inspect now three possibles cases for p. 



If < p < i(l 



2v- p 
which in turn implies that 



then 1 



pp 



> i so that S;^>i A;, 



p-3+e 



p-3+e ^ 2 ^k>l ^'k < +°°' 



L' 



^^<Cp,p,T[logN' 



k=i 



3-E 
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Furthermore 1 — § > ^—^ . 

• If i(l - ^) < p < i, then, 1 - £ < 3^ and 1 - = | so that 
Ek>i\ <+oc 

• If i(l - ^) = p < i, then 1 - £ = 3^ and 1 

w 1 ££_ 

Ek=iK ""'^^ < Cp,p,rloglog7V' (keep in mind L' 
Hence, for every e' £ (0,e), 

X:^ = o((iogiv'r^). 

fc=l AT/ 
As conclusion, we get that 



which completes the proof since log(l + A^') > ^ log (which implies log N' > 
ilog(Ar/2)). ^ 

Proof of Theorem 14.41 (a) Owing to the monotonicity of the L^-norms, 
it is enough to show that, the announced convergence holds for every q > 2 
and every p > or equivalently for every p > 2 and every q > . This 
statement follows for the ^-Holder (semi-)norm follows from Theorem l4.3r c). 
Indeed W satisfies the Kolmogorov Kp^g with 9 = p/2 — 1. On the other 
hand, it follows from [13| that, for any sequence of (Voronoi) quantizations 
at level N converging in L2 (P) toward W, this convergence also holds 

T 

in the a.s. sense. So Criterion(c) is fulfilled. 
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3-E 
Tf\ P 



2pp -, \ 1 



X L,,(Ar')OT 



A2p-1 



irhs = 2 so that 

L,,.{N') ~ logiV'). 



< 



k 



O (log iV 



1 

2(1A^) 



(4.35) 



(b) Let q > 2. The process W — W satisfies Kp^pp^i for every p G (p) 2) 
with "Kolmogorov constants" 

C^p' = C^p,T,p,.,.(logiV)-^[(^-'')"('^)], ee (0,3). 

We wish to apply Lemma 14.21 (and the remark that follows). 

> Assume < p < Then there exists ry > such that p < p' = — rri— . 
Set 9' = ^. One cheks that ^ + | < ^ so that there exists rj' > such 
that P=^ + ^+V<^- Elementary computations show that ^ — p < ^. 
Let eE (0,3) such that ^ — P < — Consequently, Lemma (and the 
remark that follows) imply that 

pi 

and for any small enough a > 0, one my specify r/, rj' and e so that \ — p = 
^(1 — |) — a. Finally this bounds holds true for p G (0,p') since the the 
L^-norm is non-decreasing. 

> Now, if p > one checks that ^ ^ ^ ~ + 1^ • It becomes impossible 
to specify p G (0, ^) so that 0' = 2<0 = ^p — 1 and 1 — p > So the 
same specifications as above lead to 

< Cg,^,r,',e,T(logiV)"^ 

v' 

which yields the announced result. <> 

5 Convergence of stationary quantizations of the 
Brownian motion for the pg-Holder distance. 

In view of what will be needed to apply this theorem to the Brownian 
motion and its functional quantizations, we need to prove a counterpart 
of Lemmas 14.21 and 14.31 for W^^. However, for the sake of simplicity, by 
contrast with the previous section, we will only deal with the case of the 
Brownian motion and its stationary quantizations. 

The main result of this section is the following Theorem. 

Theorem 5.5. Let q > 2. 



W-W^\\g,H0l 



W-W\\g^H0l 
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(a) Let {W^)n>i be a sequence of stationary quadratic functional quantizers 
of a standard d- dimensional Brownian motion W defined by \2.11\) or (2.11) 
converging to W in a (purely) quadratic sense, namely \\ \ W — W^\i2 ||2 — >• 
as N ^ oo. Then, 

yq>2, Vp>0, /9g(W,W 



as N ^ oo. 



(6) Let q > 2. Assume that, for every N > 1, is an optimal product 
quantization at level N ofW. Then, for every q > 2 and every p > 0, 



IW^ - W 



2,Af I 



If, Hoi 



= o((logAr)-t""'°(^(^-|)'p)+"), Va>0, 
so that, finally, 

;0,(W,W^) =o((logiV)-|--(7(i-f)'i)+-V Va>0. 

(c) If r > I , then 



Py(W,wL"'^"j) = o(n 



Va > 0, P-a.s. 



Note that the result of interest for our purpose (convergence on multi- 
dimensional stochastic integrals) corresponds to (2,3). The proposition 
below appears as the counterpart of Lemma 14.21 on the way to the proof. 

Proposition 5.3. Let p > 2. 



Let be defined by ( |^.i^7| ). For every 9' G {0,p — 1), there exists 



(2) 

random variable £ L^ such that 



"-a.s. Vs, te[o,r], 



|<4f)|t-s|-. 



(6) Let 



f\wi-Wl)dWi] (u;), s,te%T],s<t, 

Js / i,j=0,...,d 



where W = is a stationary quantization ofW( the integration holds in 
the Stieltjes sense). Then, for every p > 2 and every 9' £ {0,p — 1), 



'-a.s. Vs,tG[o,r], \wff\<E{zi;i,\g^)\t-s\ 



A2) 



34 



(c) Let W;f = Wlt- W;f where W = is now an optimal quadratic 
product quantization of W at level N. Then, if p > ^, for every 6' G 
{0,p{p + |) — 2), for every e G (0,3) and every 5 > 0, there exists a real 
constant Cp^p^T,d,6,S > such that 



|W 



sup 



2, TV I 
s,t I 



<Cp,p,T,d,eA^OgN) 



:,t6[0,T] \t-s\^ 

Proof, (a) The random variable Z^^ of interest is defined by 

4^=|E2"'-p,|w|,i. 

Let s,tG [0, T], s < t < s + We know from the multiplicative tensor 
property that, for every uG [s,t], 

and that, for every i, j G {0, . . . ,d}, 

\wiu^wi,\<^{\wij' +\wi/). 

2 

To evaluate sup^^j^ ^^jijIW^'J, we may restrict to dyadic numbers owing to 
the continuity in {s,t) of W^^. As a consequence, we have, still following 
the classical scheme of Kolmogorov criterion 



sup IW^J < max |W' 



m>0 



''k 



+ max \W^n 



Now 



and 



0<fc<2"+'"-l ' ""fe ' fc+l 



E max \W%+rr, .n+mf < E | W?„+^ |^ 

0<fc<2"+'»-l '•fc ''fc+l 



2m+n ^ 

E max I W. 

0<fc<2"+'" — 1 fe 'fc+l ^-^ '^+1 
^=0 
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where the norms | . | are the canonical Euchdean norms on the spaces A4((d+ 
1), {d + 1)) and respectively. 

It is clear that, for every i ^ j , i, j > 1 and every t > s, 

'V^ - wi)dwi 



I Hp 



< 



< Cp\t'-t\ 



(w: - wiYdu 



whereas 



\\\Wt'-Wt\\ = \t' 



Noting that PFt" = t and, if i = j, 1 < i < d, W^'f = \{ 



'illlp = Cp4\t' - t\. 

-iwi - wif 

shows that the above upper-bound still holds for i = j and i or j = 0. 
Consequently, we also have 

\\wlfl<Cp,,\t'-t\. 

Consequently 



2^71 -\-m -j^ 



E max \W% + m ,n + mF < Cpd > 



so that 

.(2) 



k=0 



2"- 



(1 ,^9("+™)(l-p) 



-oo 



n>0 m>0 



n>0 



since 9' < p — 1. 

On the other hand, one has obviously 



sup 



^ (2) ^ ^^^^ 



s,te[0,T],s^t 1^ _ " 

Lemma l^^ g) applied to W (which satisfies {Kp E_i)) yields for every 6' £ 
(0, 1 - 1) the existence of Z^^) € LP(P) such that 



sup ^ < Z^, 

:,te[o,T],s<t |t _ 



a.s. 
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As a consequence, combining these two results shows that, for every 



^ 2p 



where Z^^) is related to 6*' = 2 g (0, | - 1) and Z^^) is related to 6^' = ^ G 
(0,p-2). ^ 
(6) If i 7^ J, < i,j < d, it follows from Proposition 12.21 that W^'J-''^ = 
IE(W^;f^ I ^Jv) where Gn = (7{W) and T?^ = {W'^^)i<i<d is an optimal 
product quantization at level N (which means that for each component W^, 
VF*'^ is an optimal product quantization at level A^' = [A^^J). 

When i = j>l, |W^;f'^l < ((^7 - Wj)^ | Qn)). One derives that 

< ¥.\^^\gA<^({zf})'-\QN\. 



t — S P \ t — S h' 



When i = j = 0, W^'**'^ = W^'** = i(t - s)^ 
(c) In this claim, the random variable Z^?^'^ of interest is defined by 



■^e' = ^2^2 f sup |W 

n>0 s<<<s+^ 



and we aim at showing that it lies in LP(P) with a control on its L^-norm 
as a function of 
when s < u < t 



" 2 A'' 

as a function of A^. One first derives for W^'^ the straightforward identity 



where 



wff = w^;.^ + wjf + iyil. 



= {Ws,u-W,'^j0Wu,t + W,'^,^{Wu,t-W^^t) (5.36) 
with Wr^s '■= Wr — Ws if r > s, etc. One derives from ()5.36p that 
|W I < 2 > max |W 

0<fc<2"+™-l *ft '*fc+i 

m>0 



+ 2 ^ ^ max I M^,„ + m ,n+m ~ ^,n+m ,n + m 11^/"+™ +" + "11 (5-38) 

0<fc<2" + '"-l *fe 'fc '*fe + l fc' ' fc' + i 
m,m'>0 , / 

~ 0<A:'<2"+'" -1 

+ 2 \ ^ max |I^,„'+m ,n+m ~ ^,n'+m ,n + m I I | • (5.39) 

■"Hl^ 0<A;<2"+'"-l *fc ''fc+1 *fc '*fe+i fc' ' ft'+i 



m,m'>0 , I , 

- 0<fc'<2"+™ -1 
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where we used that |n t)| < 

We will first deal with deal with the first term in ()5.37p . We note that 



E 



0<fc<2 



max IW^'^ \P < V EIW^'^ 

lildA I „ + m ,n + m I — / , ^| ,n + m ,n + n 



0<fc<2"+™-l 



Let s, t G [0, T], s < t and i, j G {1, . . . , d}, i 7^ j. One checks that the 
following decomposition holds 

ct rt 



w 



s.t 



{A) (B) 

Let us focus on (A). First not that, owing to Proposition 12.11 applied 
with = a(W^, u G [0,r], Wi, s < t), 

ft 



n>l 



U 



du. 



Using that and are independent, one derives that [A) is the terminal 
value of a martingale with respect to the filtration o"(^^, k < n, Wl,0 < u < 
T), n > 1 so that combining B.D.G. and Minkowski inequalities yields, with 
the notations of Lemma 14.31 



\n>l ^-^^ 



cos 



u 



du 



< a 



BDG 



,n>l 



Wl_„cos 



'Xr, 



du 



p 

2\ 2 



2\ 2 



where £,n = £,n — £,n" and A^i, . . . , Nn, ■ ■ ■ denote the optimal bit allocation 
of an optimal quadratic product quantization at level N' (keep in mind that 
Ni^ = 1, k > Lg{N') and Ni - ■ ■ Ni^^f^^qi-^ < N' [B scalar Brownian motion). 
Now an elementary integration by parts yields 



Wl .„ cos 



u 



du 



sm 



t 



sm 



dWl 



so that, for every (0, i), one checks that, owing to the BDG Inequality, 



Win COS 



du 



— Cp^pXn^ \t 



S\ 2 



+P 
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Finally, for every eG (0,3) 

mil <Cp,T,p J Y.^n''\\^n\\l\ \t-s\^ 



, ra>l 



One shows likewise the same inequality for (B) once noted that 



i: 



wi:^ cos 



du = E 



Wlu cos 



u 



which implies 



/' 

J s 



n'^cos 



du 



< 



Win COS 



u 



du 



Consequently, for every e£ (0,3), 



|W 



s,t 



+p 



n>l 

< Cp^T,p,d,e{^OgNy 



If i = j > 1, then 



W 



2,ii,N 
s,t 



N^2 



SO that, using again Holder Inequality, 



(5.40) 



and one gets the same bounds as in the case i ^ j. 

If i or j = 0, one gets similar bounds: we leave the details to the reader. 
Finally, one gets that, for every i, j £ {0, . . . , d}, 



W^fll, <Cp,p,T,d,.,5(logiV) 



-{|-P)A 



2(p+S) U _ s\2+P, 



By standard computations similar to those detailed in Lemma 14.21 we get 

Vll max IW^'I^ \\l<Cp,,,TAeA^°sN)~^^^ 



m>0 
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Let us pass now to the two other sums. We will focus on (|5.38p since 
both behave and can be treated similarly. 



max |^^,n + m,7i+m ~ ^.n'+m ,n + m W^/-' 

0<fc<2" + '"-l 'fc + 1 'fc '*fe + l fc' ' fc' + l 

0<fe'<2"+'"'-l 

— l^,n'+m ,n + m ~W ,ri + m ,n + m \^\W_j.n + m ^n + m \^. 
' '-fe ''^fc + l ''fe + 1 fc' ' fe' + l 

0<fc<2"+'"-l 
0<fc'<2"+'"'-l 

Now for every s,u,t€ [0, T], s < u < t, it follows from Holder Inequality 
that 

\\\Ws,u-W,'!u\\Wu,t\l < \\Ws,u-W,'!u\U\\WuAl,^^,,s, 

< c,4Ws,u-w,'!u\l+,\t-u\l 

Using Inequality (j4.33p from Lemma 14.31 we get for every p > 2, every 
/9G (0,^), every eG (0,3), and every s, tG [0, T], s < t, 



Now, 



E max \W '+rn ,n+m ~ W '+m ,n + rn\^\W,n+rn ,n + m\^ 
0<fc<2" + ™-l 'fc ''fc+1 'fc '*fe + l fc' ' fc' + i 

0<A:'<2"+™'-l 

and we use that p > - and p > 2 to show that 



\^ max II I ^^,n+m ,7i+m ~ ^.n'+m ,n + m 1^1 ^^+" + "1 I |L 

■Hl^ 0<fc<2"+™-l 'fc ''fc+1 'fc '*fe+l fe' ' fc'+l 

m,m'>0 ~ ~ , / 

~ 0<fe'<2"+™ -1 

Finally, we get 

E{zf'y < c,,,,^,,,,,,(iogiv)-(i-^)-fe 

as soon as 0' G (0,0) with = p(/3 + i) — 2. Now, it follows by standard 
arguments that 

sup |w I < zy \t-s\p 
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so that, finally 



|W 



sup — 

s,te[o,T] u 



2,N, 
s,t I 

I P 



2(p+i) 



Now, we are in position to prove the main result of this section. 

Proof of Theorem 15.51 (a) Given Theorem 14.41 this amounts to proving 
that ||W^ — W^'^ll? jj^i converges to in every LP(F). This easily follows 
from Proposition 15. 3f a)-(6). 

(b) We inspect successively four cases to maximize min(l — p, ^) in p when 
it is possible. 

> qe{2, 4) and p < 2^^- Let p' be defined by ^ = ^^f^ + f with a > 
small enough so that p' > p and h + 7. < \- Then set p' = 7 + 4 — i + f 

l} (J ^ ' (j l) ^ ^ 

(note that p' > ^). One checks that i - p' = 1 - 2(^ + J) = f (1 - |) - a G 
(0, 2u7ls) ^ 5) least for any small enough a, 6 = d{a,q) > and e = 



2(p'+<5) " 2 

e(a,g) > 0. Now, Proposition 15.3( c) applied with 0' = ^ ^ f vr 1 2 
yields the announced asymptotic rate for 
since LP(P)-norms are non-decreasing in p. 

[> q £ (2,4) and p > nf'^lo) • ^^^^ same specifications as above for 



2 ' 



, p < p', 



p but with p' = p. Then 1/2 — P > ^ and choose e = e{q,a) > and 
(5 = 6{q, a) > small enough so that 2{p+5) — ^ "I" 

> q £ [4,20/3). Then 2(1-2) ^ checks that the cases p £ 

(2, 2(1-2) ) ^'^'^ ^ ^ [ oj-J^oi ' solved as above. If p > (hence 



> 5), no optimization in p is possible i.e. any admissible p satisfies ^— p > ^. 

> g> 20/3 i.e. > ^. If p < set such that ^ = ^ + aV2, 



2(q-2) ^ g-4- -^'^ ^ g-4' ^ ''"^^^ p' 2g 

2 I ^ _ 1 I a 
q ~^ p' 2 2 ■ 



a' > small enough and p' = | + — ^ + f- Doing as above yields 
min(l — p, ^) = I + a for an arbitrary small a > 0. Note that this quantity 
is greater than |(1 — |) + a (so in that case our exponent is not optimal). 
If P > we proceed to no optimization in p. 

(c) This is a consequence of Borel-Cantelli's Lemma by considering p > 
Now we conclude by proving Theorem [ 
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Proof of Theorem l3.2l First we check using Proposition[0]that pg{W^ ,0) 
and pq(W,0) are a.s. finite since they are integrable. Now we may apply 
Theorem 13. II which yields the announced result. ^ 
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about several versions of this work and F. Delarue and S. Menozzi for initiating our 
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Appendix: Functional conditional expectation 

Let (y/)/e[o,T] be a bi- measurable process defined on a probability space {Q.,A,'¥) 
such that ^ 

/ E(r/)dt < +00. 
Jo 

One can consider y as a random variable Y : {Q,,A,f) := L^{[0,T],dt) and 

more precisely as an element of the Hilbert space 

{ft, A,¥) := ^^Y : {fl,A,¥) ^ Ll, E\Y\l2^ < +00} 

where |/|^2 = Jq .f^{t)dt. For the sake of simplicity, one denotes ||y||2 := y^E jF II2 ■ 
If ;B denotes a sub-a- field of A (containing all P- negligible sets of A) then L^a B, P) 
is a closed sub-space of A, P) and one can define the functional conditional 

expectation of Y by 

E(y|B) := Proji2^(^,B,r)(y). 

Functional conditional expectation can be extended to bi-measurable processes Y 

such that liyjli := Ejyjj^i^ < +00 following the approach used for R'^-valucd 
random vectors. Then, E(y | S) is characterized by: for every S([0,T]) (g) B-hi- 
measurable process Z = (.^t)te[o,T]) bounded by 1, 

E / ZtYtdt = E f Zt E(y I B)t dt. 
Jo J a 

In particular, owing to the Fubini theorem, this implies that as soon as the process 
(E(yt I B))tg[o,T] bas a B([0, T])(g)B bi-measurable version, the functional conditional 
expectation could be defined by setting 

E(y I B)t{w) = E{Yt I B)(w), (w, t) e O X [0, T]. 
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Examples: (a) Let B := <T{J\f^, Bi, i G I) where {Bi)i^j is a finite measurable 
partition of O such that P(-Bi) > 0, i£ I. 

(b) Let Y := {Wt)te[o,T] a standard Brownian motion in R"^ and let B := a{Wti , Wt„) 
where = to < h < . . . < tn = T. Then 

Vt e [tfe,tfe+i), E{W\B)t = Wt, + - mj. 
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